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ABSTRACT 


A new analytical model for chemoviscosity variation during cure of 
thermosetting resins has been developed. This model is derived by modifying the 
widely used WLF (Williams-Landel -Ferry) Theory in polymer rheology. Major 
assumptions involved are that the rate of reaction is diffusion controlled and is 
linearly inversely proportional to the viscosity of the medium over the entire 
cure cycle. The resultant first order nonlinear differential equation is solved 
numerically, and the model predictions compare favorably with experimental data 
of EPON 828/Agent U obtained on a Rheometrics System 4 Rheometer. It has been 
shown that the model can describe chemoviscosity up to a range of six orders of 
magnitude under isothermal curing conditions. The extremely non-linear 
chemoviscosity profile for a dynamic heating cure cycle can be predicted as 
well. The model is also shown to predict changes of glass transition temperature 
for the thermosetting resin during cure. The physical significance of this 
prediction is unclear at the present time, however, and further research is 
required. From the chemoviscosity simulation point of view, the technique of 
establishing an analytical model as described here can easily be applied to any 
thermosetting resin. The model thus obtained can be used in real-time process 
controls for fabricating composite materials. 
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CHEMOVISCOSITY MODELING FOR THERMOSETTING RESINS - I 


I. INTRODUCTION 

Perhaps the single most important property of a polymer with regard to 
specifying its processing characteristics is its viscosity, which governs the 
resin flow characteristics. For thermoplastic materials, the viscosity is 
influenced by local flow geometry and can vary with temperature and shear rate. 
Viscosity control becomes more critical and difficult in the processing of 
thermosetting resins, because of the onset of chemical reactions and the 
generation of heat during curing which causes the viscosity to vary with time. 
The term chemoviscosity refers to the variation of viscosity due to polymer 
chemical reaction. The study of chemoviscosity is generally called 
chemorheology . 1 

In a typical autoclave operation for the fabrication of composite 
materials , 1 the viscosity-time profile must allow first for the bubbling off of 
trapped gas, then fiber compaction with resin flow, and, finally, laminate 
consolidation under applied pressure before the resin has gelled and ceases to 
flow. At the start of a cure cycle, polymerization begins and the polymer chains 
lengthen. The increase in viscosity of the resin due to polymerization is, 
however, largely offset by the increase in temperature which is introduced 
through the cure cycle and the heat of reaction. Consequently, a drop in resin 
viscosity of two to three orders of magnitude at the initial stage of cure is not 
uncommon. As the resin continues to polymerize, crosslinks are formed and the 
viscosity of the system starts to increase at a faster rate. Finally, the rate 
of increase of viscosity approaches infinity at gelation. Any chemoviscosity 
profile is therefore closely related to the reaction kinetics of the resin system 
and the cure cycle (temperature profile) during processing. 



Numerous researchers have devoted considerable efforts to establish an 
analytical model for the chemoviscosity growth profile of thermosetting resins. 
The most common approach to this problem has been by empirically specifying the 
viscosity linearly as an exponential function of reaction time 2 " 5 or, with the 
aid of DSC thermal analysis, the degree of cure . 6 " 7 The uses of such relatively 
simple models were apparently motivated by the empirically observed nearly linear 
viscosity growth profiles under isothermal cure conditions. A correlation 
between such isothermal cure model parameters and the curing temperatures can 
then be established and used to simulate the chemoviscosity growth in a 
non-isothermal curing condition such as that encountered inside the autoclave. 

As expected, the nonlinearity of the viscosity-time profile associated with the 
advancement of resin cure cannot be accounted for satisfactorily by such models. 

The objective of this report is to describe the applicability of the well 
established Williams-Landel -Ferry (WLF) theory in polymer rheology to the 
chemoviscosity modeling of thermosetting resins. The viscosity -time profiles 
under both isothermal and dynamic heating conditions are discussed. 
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LIST OF SYMBOLS 


c 


C 


i.C 


2 


Cl 


k 0 .n 0 

k l»2,3,4 

*T 

K 


t 


tf 

T 

To 

T g (t) 


) 


Heating rate (°K/min) 

Constants in WLF Theory, Eq. (1) 

Initial value of c\ selected for the recursion formula 
Eq. (6a) 

Constants for the parameters in Eq. (6b) 

Parameters defined in fourth-order Runge-Kutta method 

Rate constant of reaction at temperature T 

= (dT g /dT) T 
a 'g 

Time in minutes 

Time (min) for complete cure of resin 
Temperature (°K) in cure cycle of thermosetting resins 
Starting Temperature (°K) of the cure cycle 
Glass transition temperature (°K) of the thermosetting 
resin during curing 

Viscosity (poise) growth of the thermosetting resin 
during the curing temperature T 


II. THEORY 


The Williams-Landel -Ferry (WLF) theory 8 states that for 


T g <T<(T g + 100° K), 


temperature-dependent viscosity follows the expression: 


log A 

\ 


C,(T-T g ) 

C 2 + T-T g 


(1) 


where Cj and C 2 are constants. Equation (1) implies that at temperatures above 
T g all viscoelastic properties governed by the segmental relaxation rate 
(e.g., viscous flow, mechanical, and dielectric relaxations, etc.) will vary with 
temperature in essentially the same way for all polymers. The molecular 
structure effect on the viscoelastic properties of various polymers will largely 
disappear when the polymers are compared in corresponding states. The tempera- 
ture dependency of the viscosity as shown in Equation (1) has been tested and 
found valid for materials ranging from diluted polymeric systems and thermo- 
plastic melts to rubbers and elastomers. 9 The equation has also been widely used 
to describe the time-temperature superpositions for the viscoelastic properties 
of many thermorheologically simple high polymers. 10 Equation (1) can be derived 
from the semi -empi ri cal Doolittle equation which relates the viscosity to the 
free volume of the liquid. 11 As the free volume increases, the viscosity rapidly 
decreases. The two constants C x and C 2 are, therefore, related to the fractional 
free volume at T g and the thermal coefficient of expansion of the fractional 
free volume above T g . Cj and C 2 were originally thought to be universal 
constants with values of 17.44 and 51.6 respectively. It was later found that 
they are material-related. A tabulation of these material constants for various 
systems can be found in the literature. 9 
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The normal use of the WLF equation for thermoplastic materials requires that 
the glass transition temperature Tg be constant while the temperature T is 
varied for the specific polymer under study. Eq. (1) is applicable for a 
temperature T up to 100°K higher than the T g of the material. However, during 
cure of thermosetting resins, the monomers are initially polymerized and 
crosslinks are formed later. This is a system where T g (t) is changing and the 
curing temperature T is held constant, in the isothermal case. The glass 
transition temperature T g rises continuously and may eventually approach the 
curing temperature. Over the entire curing cycle, the material structure 
actually undergoes continuous phase transformations from the low molecular weight 
liquid to a high molecular weight polymeric melt, and eventually transforms to 
form crosslinking networks. If it is assumed that T g (t) of the material is 
always lower than the cure temperature T, and that (T-T g (t)) is always within 
100°K, then the WLF theory should be applicable to all of the different polymer 
structure phases during cure. Before Equation (1) can be applied to describe the 
chemoviscosity of the curing resin, however, modifications have to be made in 
order to account for the time factor. This can be accomplished by the following 
rational izations. 

It may be assumed that a segmental relaxation rate governs most rate 
processes that take place in polymers at temperatures above T g . The rate of 
reaction is often diffusion controlled, especially in the crosslinking stages. 
According to the accepted theory of diffusion controlled reactions, the rate 
constant kj is generally proportional to the diffusion constant of a reactant, 
and therefore is inversely proportional to the viscosity of the medium. 
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It is further assumed that the rate constant of reaction kj is directly 
proportional to the rate of change of the glass transition temperature Tg, 

k "l" n dT dT 
= ^T^dt^T 

Equations (1), (2), and (3) can then he combined to form a first order ordinary 
nonlinear differential equation: 


(3) 


dT Ci(T-T ) 

(dt^T = K exp tc' +T-T 9 ] ( 4 ) 

2 g 


where K = (dTg/dt)-^ is the rate of change of Tg at T = Tg, and 
is taken to be a material constant. 


With the boundary condition that as t+tf, Tg+T where tf denotes the 
time for complete cure of the resin and T denotes the cure temperature. 

Equation (4) can be solved numerically for T g as a function of curing time, 
t. The chemoviscosity of the curing resin can then be obtained by Equation (1). 

In the case of dynamic heating curing conditions, T in Equation (4) can be 
replaced by (T 0 +ct), where T 0 is the starting temperature and c is the 
heating rate in (°C/min). For a given cure cycle (i.e., temperature profile 
inside the autoclave), two material constants, n T and (dT a /dt) T , 

g y 'g 

and the two WLF constants and C 2 need to be specified before Equations (4) 
and (1) can be solved. Details for solving these equations will be discussed in 
Section V. 
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III. EXPERIMENT DESCRIPTION 


A diglycidyl ether of bisphenol A based epoxy resin EPON 828 (Shell) was 
studied. The hardener selected was a liquid amine EPON Curing Agent U, also from 
Shell. EPON 828 and Agent U were weighed separately with a microbalance which 
has a resolution of 0.001 gm. A weight ratio of 100/25 was studied. The two 
components were hand mixed for 2 minutes before being transferred to a 
Rheometrics System 4 Rheometer for chemoviscosity measurements. 

The sample was confined in the gap between two parallel plates mounted in 
the Rheometer. The top plate is motor driven about its axis while the bottom 
plate is mounted on a torque transducer for force measurement. The typical gap 
between the parallel plates was 1.2 mm. The plates and sample are enclosed in a 
heat chamber where temperature control is provided. Though the shear field in 
parallel -pi ate flow geometry is non-uniform, for the small amplitude 
displacements used, these effects can be safely neglected. Plate oscillatory 
motion was set at a frequency of 10 rad/sec, and an amplitude (strain) equivalent 
of 1 percent of the gap was used for all measurements reported. Selection of the 
strain value was to assure that the measurements were performed within the 
material's linear viscoelastic response range. 

It is well known that the viscoelastic properties of polymeric materials 
respond differently to different frequencies, and selection of the frequency 
value was made to assure that the complex viscosities measured were within the 
Newtonian region. 12 The cure mode of the Rheometer was used during testing which 
automatically programs different temperature profiles as the epoxy is cured. 

Both isothermal and dynamic heating measurements were performed. 
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IV. EXPERIMENTAL RESULTS 


Chemoviscosity increase profiles for isothermal curing at three different 
temperatures are shown in Figure 1. The first data point, measured at t = 0, 
corresponds to approximately 5 ± 0.5 minutes delay from the beginning of sample 
mixing in all cases. The initial viscosity is lower when the material is cured 
at higher temperatures. This is believed to be due to the temperature effect. 
Figure 1 shows the nonlinearity in the viscosity-time profiles, especially during 
the initial stages of cure. The profiles become more linear as cure progresses. 
From the process point of view, the initial stages of cure are far more important 
because the resin flow is easier to control. A linear chemoviscosity model in 
this region would certainly be inadequate. 

The chemoviscosity profile for a dynamic heating condition is shown in 
Figured. The temperature profile is also included in the figure. The sample is 
first heated at a rate of l°C/min from 32° C to 62° C, then the temperature remains 
constant throughout the cure cycle. The viscosity initially decreases and then 
increases exponentially. Such viscosity profiles are typical for thermosetting 
resins under dynamic heating conditions. The decrease in viscosity during the 
heat up period provides more flexibility in resin flow adjustments for the 
fabrication of composite materials in an autoclave. Precise chemoviscosity 
modeling is therefore more important under these conditions. 
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V. MODEL PREDICTIONS 


Isothermal case 

Using Equation (4), which is rewritten here, 


dT C,(T-T ) 

^dt“ W = K exp f C 2 +T-T 9 1 
2 9 


(4) 


the glass transition temperature T g (t) can be obtained, 
of the curing resin can then be calculated by 

- c i( T - T 0 ) 

n T (t) = exp [- +T _ T 9 ] 

g 2 g 


= (n T • K)/(dT /dt) 

g 9 


The chemoviscosity 


(5) 


The two material constants K and n T have been assigned values of 

9 

2 x 10- °K/min and 1 x 10 8 poise, respectively, for all the isothermal model 
predictions described below. Physically, these two material constants represent 
the material properties of the resin at the completely cured state. 

Constant C x in the WLF theory is assumed as a function of T g by the 
following recursion formula 


with 


c i = ^(t* /T; +1 ) n 


TgVt) = T g (i At) , i = 0, 1, 2, 


(6a) 
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and 


k o t 

n = n e 
o 


(6b) 


while C 2 is assigned the "universal constant" of 51.6°K. The change of as 
the resin curing progresses (Eq. (6a)) is required by noting that values of 
vary for different polymeric systems ranging from dilute solutions to cross 
linked rubbers. 9 

The initial condition, T g (0) at t = 0, is difficult to specify but is 
required in order to solve the differential Equation (4). Theoretically, 

Tg(tf) should be equal to the cure temperature T when the time tf for 
the completely cured stage is reached. However, in practice tf may well 
approach infinity. In the present study, T g (0) is arbitrarily chosen, and the 
simulated results for T g (t) that were obtained can only describe the relative 
changes of T g as a function of time for the resin under study. 

The constants used to solve Equations (4), (5), and (6a, b) for the 
isothermal curing of the EPON 828/Agent U system with weight ratio of 100/25 
reported here are summarized in Table 1. 


TABLE 1. Constants used in the theory for isothermal case. 
(EPON 828/Agent U, 100/25 weight ratio) 


Cure Temp. 
T (° K) 

Tg(0) 

(°K) 

K 

(°c/min) 

% 

(poise) 

C° 

c 2 

n o 

k o 

305 

280 

2*10- ? 

1-10* 

i 

41.5 

51.6 

1 

0.165 

310 

280 

1 

O 

r- H 
• 

CM 

CO 

G 

i— f 

• 

r— 1 

38 

51.6 

1 

0.245 

317 

280 

2.10-A 

1.10 8 

37 

51.6 

1 

0.385 
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Equation (4) is solved numerically by a fourth-order Runge-Kutta method 13 as 
modified by Gill: 


T g 1 + I * V ["> + 2 <1 - b M + 2 (1 + b K 3 ♦ M 

/ 2 ic 

k i ■ fit 1 , Tg 1 ) 


k 2 = f(t'+£, Tg 1 + £k,) 

k 3 = f(t* + £■, T * + (- i + ^)hk x + (1 - i)hk 2 ) 

y n n 

k 4 = f ( t ^ + h, T "* - — hk 2 + (1 + — )hk,) 

9 ft n 


where h is the time increment and 


(7) 


c i( T - T J 

f(t, T ) = K exp [£_*-] (8) 

2 g 

This algorithm was programmed in Basic and solved on a Hewlett-Packard HP-9836 
computer. The program listing is included in the Appendix for reference. 
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Numerical solutions are plotted in Figure 1 with the experimental data for 
isothermal curing at three temperatures. Good agreement between the data and 
theory for all three cases are noted. The capability of the model to describe up 
to six orders of magnitude in chemoviscosity range is significant. Numerical 
computations were performed for different time increment (At and h in Eq. (6a) 
and (7), respectively). It was found that the simulated chemoviscosities based 
on At = 0.5 and 0.25 minutes, respectively, differ by less than 1% for all three 
temperatures. The simulated initial viscosities for various temperatures are the 
result of different assigned values of C? and n 0 , while the simulated viscosity 
increase rates are dictated by the assigned k 0 . Variations of the computed 
are plotted in Figure 3. It is noted in the literature 9 that Cj varies from 11.2 
for polybutadiene rubber to 34.0 for Methyl (atactic) polymer. The variations of 
C x shown dictate structure changes of the material under investigation, and are 
approximately equivalent to a variation from 0.011 to 0.087 for fg, the 
fractional free volume at Tg, of the curing resin. 

Calculated values for Tg(t) for the curing resin are also shown in Figure 
3. It can be noted that changes of Tg are faster for higher cure temperatures 
during early curing stages. Higher cure temperatures yield higher Tg. After 
the initial curing period of approximately 20 minutes, values of Tg quickly 
level off; the chemoviscosities n(t), however, keep increasing as can be seen 
from Figure 1. Similar situations exist in the relationship between n(t) and 
a(t), the degree of cure determined by DSC thermal analysis, as reported 
experimentally by Lee et al. 7 These observations imply that determinations of 
gelation points of the thermosetting resins by viscosity data alone are 
questionable. Further research in this area is planned. T g (t) in Figure 3 
increases very slowly after the initial 20 minutes, and may take a very long time 
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to reach cure temperatures of 305, 310, and 317°K, respectively. This is in 
agreement with the fact that low curing temperatures are employed. 

The analytical model developed here (i.e., Eqs. (4), (5), and (6a, b)) 
possesses three material constants and four adjustable parameters. In describing 
the chemoviscosity under isothermal curing, however, only two parameters (C? and 
k 0 in Table 1) need to be adjusted in order to account for the differences in 
cure temperatures. It can be noted from Figure 1 that the model predictions are 
in better agreement with the experimental data than the empirical 4-parameter 
model used in the literature 3 . The values of Cf and k 0 used follow the 
Arrhenius relationship shown in Figure 4. The Arrhenius correlation makes 
chemoviscosity predictions possible for this material at temperatures outside the 
range studied here. 


Dynamic Heating Case 


In this case cure temperature changes with time. The value of T in 
Equations (4) and (5) is replaced by (T 0 + ct), where T 0 is the start temperature 
and c is the rate of increase of the temperature. Equation (4) becomes 


dT 


Ci(T +ct-TJ 

dt 1 = K exp [r 2 vr- + ct--T ] 

1 o g 


(9) 


When Equation (9) is directly applied to the cure cycle shown in Figure 2 
using the material constants listed in Table 1, the predictions as plotted by the 
dashed line are relatively poor. The drop in viscosity in the initial cure 
portion is overestimated, and the rate of viscosity increase as the cure 
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progresses is not as high as those experimentally measured. By selecting the 
following new set of constants: K = 9.2 x 10~ 7 , C° = 41.0 and 

n Q =1, k Q = 0, for 0 < t < 15 mins 

= 0.00012, = 0.62, 15 < t < 21 

= 0.00012, = 0.59, 21 < t 

the simulated results shown by solid line in Figure 2 are in very close agreement 
with the experimental data. 

When numerical computations were performed for two different time increments 
(At = 0.5 and 0.25), the values of calculated chemoviscosity differed by less 
than 1%. The variations of Tg(t) and Cj(t) during cure are shown in Figure 5. 

The behaviors of Tg and Cj are very similar to those discussed before under 
isothermal cure cases. 

The theory was also compared with dynamic heating cure data reported by 
Tajima and Crozier. 1 The material used was composed of two epoxy resins, 
TGMDA-MY720 of 69 wt. percent and SU-8 of 9 wt. percent, and an amine hardener 
0DS of 22 wt. percent. The cure cycle for this resin together with the 
chemoviscosity data are reproduced in Figure 6. Material constants selected for 
the model are as follows: 
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T (0) = 284° K, K = 2 . 10" 10 °k/min 

n T = 1 x 10 10 poises, C° = 31.6, C 2 = 51.6°K 

9 


n 


o 


0 . 1 , 


k = .001, for 


= 0 . 0002 , = 0 . 1 , 

= 4 x 10- 5 , = 0.1, 

= 1 x 10- 5 , = 0.115, 


0 < t < 71 min 

71 < t < 111 

111 < t < 121 

t > 121 


Excellent agreement was demonstrated between the model predictions and 
experimental data. Variations of Tg(t) and Cj(t) during cure are also plotted 
in Figure 7. The behaviors of these quantities are similar to those discussed 
before for the Epon 828/Agent U system. 


Equations (6a) and (6b) have been incorporated in Equation (4) or (9) to 
simulate the chemoviscosity of thermosetting resins under various cure cycles 
discussed so far. Equations (6a) and (6b) are, however, purely empirical. Other 
forms can also be successfully used. For example, consider the following new 
equation: 


T. k o 


C l = C 1 ] [1 + logl-f) ] 


and the recursion formula for k 0 as 


i i + i i i it i + l ,-r i \ 

k n = x (T /T ) x n A 
o o g g ' o 


(10a) 


(10b) 
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Material constants selected for this model are as follows: 


T (o) = 46° K, K = 5 x 10” 14 °K/min 

n T = 1 . 10 13 poises, C° = 32.44, C 2 = 51.6°K 
9 

while 

n Q = 0.98, k° = 0.19, for 0 < t < 75 min 

= 1.025, = 0.19, 75 < t < 111 

= 1.008, = 0.19, t > 111 

When the temperature profile in Figure 8 is followed. Equations (9), (10a), and 
(10b) give solutions plotted as the solid curve in Figure 8. The model 
predictions again agree well with the experimental data over the entire cure 
cycle. 

Tajima and Crozier 1 used the same WLF Theory to explain their experimental 
data. However, the treatments are quite different from those described in this 
paper. They choose reference temperatures T s (t) to replace Tg(t) in Eq. 

(1). While both n Ts (t) and T g (t) are varied, the constants Cj and C 2 are 
kept constant during the entire cure cycle. The values of ny s (t) and T s (t) 
are related to the content of unreacted curing agent (DDS) where the time factor 
is introduced. Physically, continuous changes of reference viscosity and 
temperature in the WLF theory imply continuous changes in molecular structures of 
the material. This accounts for the good agreement between their model and the 
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experimental data. The approach reported in this paper involves only one 
experiment (i.e. Rheometry) and the model is much easier to establish. 


Mussatti and Macosko 14 reported steady shear viscosity data on EPON 828 
catalyzed by 1% by weight of triethanolamine (TEA) under isothermal curing at 
100°C. Viscosity measurements were made with cone and plate flow geometry of a 
Rheometrics Mechanical Spectrometer. Rotation was started and stopped 
periodically as the reaction proceeded, and the resulting torque was recorded 
until the sample broke up. The viscosity data obtained in this experiment are 
reproduced in Figure 9. Attempts have also been made to fit the data with the 
model, i.e., Eqs. (10a, b) for the parameters C L and k 0 have been used together 
with the model Eq. (4). Material constants and parameters selected for the model 
are as follows: 


and 


T (0) = 150° K, K = 5.0 x 10" 14 °K/min 


n Tg = 

1 x 10 

13 poises, Cj = 38.44, 

c 2 = 51. 6° K 

bo = 

1.01 , 

k 0 = 0.1, for 0 < t < 

50 mins 

= 

1.045, 

= 0.1 

50 < t < 

60 

= 

1.022, 

= 0.1 

t > 

60 


It can be seen in Figure 9 that the theory can also describe the steady shear 
viscosity data very well. Variations of Tg(t) and C L (t) during cure are 
plotted in Figure 10 for reference. The behaviors of these quantities are again 
similar to those discussed before for other resin systems. 
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VI. CONCLUSIONS 

It has been shown that the WLF theory, Equations (4) and (9), with the 
appropriately selected material constants can accurately predict the 
chemoviscosity of thermosetting resins. The success in applying the theory to 
an entire cure cycle for both isothermal and dynamic heating conditions indicate 
valid assumptions for chemoviscosity model development. It is surprising to note 
that the simple assumption of Eq. (2) (where viscosity is assumed to be linearly 
inversely proportional to the rate constant of reaction) incorporated with WLF 
theory can account for the chemoviscosity increase over an entire cure cycle. 

The demonstrated capability to describe up to six orders of magnitude of change 
in the chemoviscosity range with the model is significant. The model can also 
predict the change of glass transition temperature T g ( t ) during cure. 

The physical significance of the predicted T g (t) which results from the 
assumption of Eq. (3) in the theory is unclear at the present time. A thorough 
understanding of the kinetics of chemoviscosity, the degree of cure, the 
stoichiometry and mechanism of polymerization can only be achieved by a combined 
use of several analytical methods, such as Rheometry, DSC, TBA, and 
spectroscopy. Such tasks remain to be resolved in the future. 

The model as described has considerable value in autoclave process 
monitoring applications. An analytical model with the appropriate parameters 
determined from a single set of experimental data (e.g., via Rheometry) can 
easily be established for any thermosetting resin. The established model, which 
has the capability of describing the chemoviscosity accurately, can be used in 
real-time control of automatic composite processing equipment when an appropriate 
sensing device can be identified. 
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APPENDIX 


Shown herein is a listing for solving Eqs. (4), (5), (6a, b) by means of 
a Runge-Kutta algorithm. Hewlett-Packard Basic Extension 2.0 language is used. 
The program is run on a HP-9836 computer. 
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10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 


! 2/8/84. “DIFFER? AO” 

! 2/15/84. REVISED"'! TO ACCOMMODATE BOTH CONST AND NONISO-THERMAL CASES 

i 

• THIS PROGRAM FIXES FINAL CURING TEMP » T AT TIME-TFINAL 

! ADJUSTABLE VARIABLES ARE K.A.B FOR THE NON- ISOTHERMAL CASE 

• NUMERICAL SOLUTIONS ARE PROVIDED BY SHOOTING ME T HOD IN- 

! CORPORATED WITH THE RUNGE-KUTTA 4TH ORDER SCHEME TO SOLVE 
! THE DIFFERENTIAL EQUATION 

• IT HAS BEEN ASSUMED THAT K-<DTG/DT> AT TG IS A CONSTANT; 

! THAT A<I»-A0(I)«<TG(I>/TG<I+1 >KN, A(I>-C1 

! 8-51 .6. B-C2 IN WLF THEORY 

! CONST AND NON-ISOTHERMAL CASE 

t » 

COM T imechgO . Timechgl . Tinechg2. T imechg3 .Cc . AO 

COM TgOtemp.TgOiast.TgOnou.Alast .Anow 

COM Nn.N0.N1 .N2.N3,Kn1 .Kn2.Kn3.Kn4,KnS.Kn6.Kn7.Kn8 

i 

TgO-284 
K-2.0E-1 0 
AO-31 .6 
N0-.1 
Kn2- .001 
T imechg0-71 

» 

N1-.0002 
Kn4- . 1 

T imechgl - 1 1 1 

* 

N2-. 00004 
Kn6- . 1 

Timechg2- 121 

i 

N3-. 00001 
Kn8- .115 
T imechg3-2000 

'let a-ao 

LET B-51.6 
Cl st- 1 . 

C2nd-.8 
Time 1-70 
T ime2-1 1 5 
Tstart-338 
Temp-Tstart 
Tf inal-1000 
Xxtg=1 .E+1 0 

i 

H-.5 

* 

'PRINT ”DIFFER3_AO_CROZIFR2“ 

PRINT 

PRINT “SOLUTIONS FOR TSTART(K' ; Tstart :“ TFINAL<MIN)-“:Tf inal 

PRINT “A0-” : AO : ”B-“ ;B : “K-“ ; K ; “H-“ :H 

PRINT “TEMP PROFILE IS C1ST-" ;C1st ;“ AT TTME1 -" ;T imel 

PRINT " AND C2ND-”;C2nd:“ AT TIME2-" : Time2 

PRINT “TIMECHG0-” : T imechgO ; “NO-” : NO ;“FN 1 ; Knl : “KN2-" ;Kn2 

PRINT “TIMECHG1 •“ ; T imechal ;"N1-":N1 :"KN3-":Kn3:"KN4-“:Ki>4 

PRINT "TIMECHG2-” : Timechi2;“N2-“ ;N2 ; "KN5-" :<n5 : "KN6-" ;Kn6 

PRINT "TIMECHG3-" ; Timechg3 ;“N3-“ ;N3 : "KN7-" ; Kn7 ; ”KN8-“ ;Kn8 

PRINT “XXTG-” .Xxtg 

PRINT 


620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
340 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1 130 
1140 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 


T imeO-O . 

Clock-1 

Dtqdt-FNValue( TimeO .TgO.A.B.TimeO.K.Tstart.Tstart) 
Xx-Xxtg*K/Dtgdt 

PRINT TiroeO.TgO." " . Dtgdt . "XX*" .Xx 

PRINT 

! TEMPERATURE PROFILE IS DEFINED AS FOLLOWS: 

IF TimeO<Time1 THEN 
Cc*C1st 

Temp-Tstart+Cc«T itneO 
Tempnou*Tstart + Cc«< T ime0 4 H> 

GOTO 880 
ELSE 

GOTO 770 
END IF 

IF TimeO<Time2 THEN 
Cc-0 

Temp-Tempnou 
Tempnoul -Terao 
GOTO 880 
ELSE 
Cc*C2nd 

Temo-Tempnou1+C2nd«< T imeO-T ime2) 

Tempnou-Tempnoul +Cc*(T imeO+H-T ime2) 

END IF 

CALL Runge_k at ta< T imeO . igO .H.Tgneu.A.B.TimeO.K.Tstart .Temp) 

TgO-Tgnew 

TgO last- TgO temp 

TgOnou-TgO 

Alast-Anou 

T imenow-TimeO+H 

Dtgdt-FNValO<Anou.B.T jpienoy , TgO . K . Temonow ) 

Xx-Xxtg*K/Dtgdt 

IF Timenow>Tf inal THEN 

GOTO 1100 

ELSE 

GOTO 1010 
END ir 

IF Timenow-Clock THEN 
Clock-Clock +5 

PRINT T imenou;Terar>nou: TgO : Dtgdt ;"XX*" :Xx :Anou 
T imeO-T imeO+H 
GOTO 700 
ELSE 

T imeO-TiraeO+H 
GOTO 700 
END IF 
END 


SUB Runge_kutta(Xi . Yi ,H . Yi 1 .A.B.Tlme.K.Tstart.Temp) 
COM T imechgO . T itnechgl , T imechg2 . T imechg3 . Cc . AO 
COM TgO temp. TgO last . ToOnow. A last . Anou 
COM Nn.NO.N1 ,N2,N3,Kn1 ,Kn2 .Kn3 ,Kn4 ,Kn5.Kn6 .Kn7 ,Kn8 
TgOtemp-Yi 

K 1 *FNValue(Xi , Yi , A .B . T ime .K . Tstar t . Temp) 

Xt-Xi+H/2. 

Yt*Yi+H/2.-K1 

Tempp-Temp+Cc»H/2. 

K2-FNValue(Xt -Yt.A.B.Time.K.Tstart.Temppl 
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1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 


Xt-Xi+H/2. 

Yt*Yi+(-.5+1 . /SQR<2 . > )«H*Ki ♦ ( 1 •-? . /$QR(2. ) >*H«K2 
K3-FNValue(Xt,Yt,A.B.Time.K.Tstart ,Tempp> 

Xt-Xi+H 

Yt- Yi- < 1 ./SQR<2.')«H«K2+n . + » . /SQR<2. > >*H«K3 
Tempp-Temp*Cc*H 

K4*FNValue<Xt . Yt .A .B.Time .K . Tstart . Tempp) 

V1-2.*<1 .-1 ./SQR<2.)> 

V2-2.-H . + 1 ./SGR<2.>> 
Yi1-Yi+H/6.MK1+V1*K2+V2*K3+K4> 

SUBEND 


DEF FNVal ue<X.Y. A. e.Time.K. Tstart. T emp> 

COM T imechgO , T ltnechgl .Timechg2 . T imechg3.Cc .AO 
COM TgOtemp . TgOlast , TgQnou . A lac- 1 . Anou 
COM Nn,N0,N1 .N2,N3,Kn1 ,Kn2 .Kn3 ,Kn4 .Kn5 .Kn6 .Kn7 .Kn8 
TO-Temp 

IF Time-0 THEN 

Nn-NO 

TgOnou-Y 

TgOlast-Y 

Alast-AO 

GOTO 1740 

ELSE 

GOTO 1490 
END IF 

IF Time<TiroechgO THEN 
Nn«N0«EXP< <Time>*Kn2> 

GOTO 1740 
ELSE 

GOTO 1550 
END IF 

IF T ime<T imeshg 1 THEN 
Nr>-N1 »EXP< T ime»Kn4 ) 

GOTO 1740 
ELSE 

GOTO 1610 
END IF 

IF T ime<Tiroecha2 THEN 
Nn-N2«EXP< T ime»Kn6 > 

GOTO 1740 
ELSE 

GOTO 1670 
END IF 

IF T ime<T imechg3 THEN 
Nn-N3*EXP(Time«Kn8) 

GOTO 1740 
ELSE 

GOTO 1730 
END IF 

t 

Anou-Alast«<TgOlast/TgOnow> ''Nn 
A-Anou 

T1-A«<T0-Y> 

T2-B+T0-Y 

T3-T1/T2 

Value-K*EXP<T3> 

RETURN Value 
FNEND 
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1840 ! 

1850 ! 

1860 DEF FNVal0(A0, BO. TimeO.TgO.K, Tempi) 
1870 A-AO 

1880 B-BO 

1890 T1-A*(Temp1-TgO) 

1900 T2 _ B+Tenip1-TgO 

1910 T3-T1/T2 

1920 Value 1 ■K«EXP<T3) 

1930 RETURN Valuel 

1940 FNEND 

1950 ! 

1960 ! 
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Figure 1. Comparisons between experimental chemoviscosity data and 

the model predictions for EP0N 828 with Agent U. Chemoviscosity 
of the resin is measured by Rheometrics System 4 Rheometer under 
three isothermal curing conditions. 
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Figure 2. Comparisons between experimental non-isothermal chemoviscosity and the 
model predictions for EPON 828 with Agent U. 

Chemoviscosity of the resin is measured by Rheometrics System 
4 Rheometer. The cure cycle (temperature profile) is also shown. 
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Figure 3. Variations of Tg(t) and C x during cure obtained from the simulation 

models in Figure 1 for EPON 828/U resin system under isothermal curing 
conditions. 
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Figure 6. Comparisons between experimental non-isothermal chemoviscosity 
data and the model predictions. Cure cycle is also included. 
The data were reported by Tajima and Crozier for TGMDA-MY720/ 
SU-8/DDS. 
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Figure 7. Variations of T g (t) and C x during cure obtained from the simulation 
model in Figure 6 for the TGMDA/SU-8/DDS resin system under dynamic 
heating conditions. 
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Figure 9. Comparisons between the experimental steady shear chemoviscosity 
data and the model predictions. Data were reported by Mussatti 
and Macosko for EPON 828/HHPA under isothermal curing at 100°C. 
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